Fractal initial conditions and natural parameter values in hybrid inflation 
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We show that the initial field values required to produce inflation in the two fields original hybrid 
model, and its supergravity F-term extension, do not suffer from any fine-tuning problem, even when 
the fields are restricted to be sub-planckian and for almost all potential parameter values. This is 
due to the existence of an initial slow-roll violating evolution which has been overlooked so far. 
Due to the attractor nature of the inflationary valley, these trajectories end up producing enough 
accelerated expansion of the universe. By numerically solving the full non-linear dynamics, we 
show that the set of such successful initial field values is connected, of dimension two and possesses 
a fractal boundary of infinite length exploring the whole field space. We then perform a Monte- 
Carlo-Markov-Chain analysis of the whole parameter space consisting of the initial field values, field 
velocities and potential parameters. We give the marginalised posterior probability distributions for 
each of these quantities such that the universe inflates long enough to solve the usual cosmological 
problems. Inflation in the original hybrid model and its supergravity version appears to be generic 
and more probable by starting outside of the inflationary valley. Finally, the implication of our 
flndings in the context of the eternal inflationary scenario are discussed. 
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I. INTRODUCTION 

The paradigm of inflation [J [2, [3, 01 is currently the 
simplest way to solve the standard cosmological prob- 
lems and explain the Cosmic Microwave Background 
(CMB) anisotropics observed so far, though other al- 
ternative mechanisms have been proposed (for a review 
see and references therein). Many models of infla- 
tion have been proposed ^6', T] , based on single field or 
multi-field potentials. If single field models are efficient 
effective models, hybrid models explore the possibility 
that the inflaton is coupled to other scalar fields, as 
first proposed by Linde _8]. When coupled to a Higgs- 
type field, inflation is realized in the so-called "inflation- 
ary valley" when the Higgs vacuum expectation value 
(vev) vanishes and the inflation end is triggered when the 
Higgs becomes tachyonic and develops and non- vanishing 
vacuum expectation value (vev). Similar models have 
rapidly been constructed in various theoretical frame- 
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works [1, [13, [m, the most popular of them being the 
supersymmetric/supergrayity versions of F-term or D- 
term inflation [l3, [ll [li, [llj . 

In the limit of sub-planckian field values, all hybrid in- 
flation models were however thought to require extremely 
fine-tuned initial field values to produce enough e-folds of 
acceleration, from the original model proposed by Linde 
to most supersymmetric versions [TtI [isj . with the 
exception of hilltop potentials which assume that infla- 
tion takes place near a maximum of the potential [l9l[20j. 
The successful initial field values were found located only 
in an extremely narrow band around the inflationary val- 
ley, or on a few scattered points away from it [2l|. This 
was considered as a fine-tuning problem for these mod- 
els since any pre-inflationary era would have to be fine- 
tuned to allow inflation to last long enough to solve the 
standard cosmological problems. This flne-tuning has re- 
cently been revisited in Ref. [l^l for the original hybrid 
model as well as for the supersymmetric "smooth" and 
"shifted" models. Using higher precision, it was shown 
that the successful initial field values are rather organised 
in intricate dense regions outside of the inflationary valley 
(see for instance Fig. 7 of Ref. ^22]). The area occupied 
by these regions was found to represent up to 15% of the 
sub-planckian field regime for the original hybrid model 
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and up to 80% for smooth hybrid inflation. The physical 
explanation of these new successful regions conies from 
the existence of an initial fast-roll phase during which the 
fields roll down the potential in a chaotic way followed 
by a climbing up of the valley and a slow-rolling phase 
back down. 

However, as discussed in Ref. ^23], these new successful 
regions appeared to depend on the shape of the poten- 
tial, and therefore on the potential parameters. One may 
wonder whether these features are a new solution of the 
fine-tuning problem, i.e. if they are robust with respect to 
the potential parameters. Moreover, Ref. [l^ did not dis- 
cuss the statistical properties of this space and the effect 
of the initial field velocities which were assumed vanish- 
ing. Finally, the popular supersymmetric extensions, F- 
or D-term hybrid models, were not studied. The purpose 
of the present paper is to quantify how these new suc- 
cessful inflationary regions are widespread in the higher 
dimensional space of all the model parameters, i.e. by 
considering not only the initial field values but also their 
initial velocities and the potential parameters. We also 
extend our analysis to the F-term hybrid model, studied 
in SUGRA 

In order to deal with a multi-dimensional parameter 
space, after having discussed the fractal nature of the 
successful inflationary regions, we introduce a probability 
measure and perform their exploration by using Monte- 
Carlo-Markov-Chains (MCMC) methods. The outcome 
of our approach is a posterior probability distribution on 
the model parameters, initial velocities and field values 
such that infiation lasts more than 60 e-folds^. As will be 
shown in the following, thanks to infiation starting "out 
of the valley" , a high number of e- folding appears to be 
generic, and favoured, in the original hybrid model for 
parameter ranges covering several orders of magnitude. 
We have also checked that such a result is not peculiar 
to a given potential by applying the same analysis to the 
more realistic two-field F-term inflation potential. This 
treatment allows us to establish natural bounds on the 
parameters (or combination of parameters) for each of 
these scenarios. 

At this point, we would like to emphasize that our aim 
is not (yet) to constrain these models with the current 
Cosmic Microwave Background (CMB) and astrophysi- 
cal data but rather to discuss in details their ability to 
generate an inflationary phase. In particular, in the small 
field limit, original hybrid model are known to generate 
a blue spectrum of scalar initial perturbations^ , which is 
disfavoured by recent CMB experiments [24] . Our use of 
this model here is motivated by its simplicity and its rep- 



^ Such probability distributions are almost independent of the cho- 
sen number of e-folds: once the field rolls down in a flat enough 
region of the potential, the total number of e-folds generated is 
always very large. 

^ This conclusion can be altered when additional couplings are 
assumed for the inflaton 1231. 



resentativity of the non-linear two-field dynamics. The 
more realistic F-term SUGRA model is in agreement with 
the current CMB data: it predicts an almost scale invari- 
ant spectrum and the generic formation of cosmic strings 
(25| . a combination which was shown to be favoured by 
observations in Ref. [2^ . 

The paper is organised as follows. In the following 
section, we discuss the fractal nature of the successful re- 
gions of inflation in the original hybrid model and define 
a probability measure over the full parameter space. In 
Sec. mil the MCMC method is introduced and we study 
step by step the effect of the initial field velocities and 
the potential parameters on the probability of obtaining 
60 e-folds of inflation. We then present the full posterior 
probability distributions of these parameters for the orig- 
inal hybrid scenario. In Sec. llVi we perform the same 
analysis of the F-SUGRA hybrid potential. Some con- 
clusions and perspectives are finally presented in the last 
section. 



II. FRACTAL INITIAL FIELD VALUES 
A. The model 

The original hybrid model of infiation was proposed in 
Refs. [iJl, its potential reads 



U(0,7^) = imV 



A 



(1) 



The field (j) is the infiaton and ■0 is the auxiliary Higgs- 
type field while A, A' are two positive coupling constants 
and m, M are the two mass parameters. Inflation is 
assumed to be realized in the false- vacuum along the val- 
ley'^ (■0) = and ends due to a tachyonic instability of ip 
when the inflaton reaches a critical value (j)c = Af -y/A/A'. 
The classical system evolves toward its true minimum 
{(j>) = 0, and (ip) = ±M whereas in a realistic scenario 
one expects the tachyonic instability to trigger a preheat- 
ing era [13, [H, [11, [Sfl [Ml, [H . 

To observe the effects of varying the free parameters in 
the dynamics of inflation, it is more convenient to rewrite 
the potential into 



02 ^2^2 



/i 



(2) 



under which Af , /i, are three mass parameters. With 
this expression, the critical point of instability now reads 



V2i^ 
M 



(3) 



Throughout the paper (.) denotes the vev of a field. 
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It is common usage to consider the effective one-field po- 
tential by restricting the field dynamics to the inflation- 
ary valley and one gets 



1+ ( - 



(4) 



B. Equations of motion 

In a flat Friedmann-Lemaitre-Robertson-Walker 
(FLRW) metric, the equations governing the two-fleld 
dynamics are the Friedmann-Lemaitre equations*, 



d Stt 



a 3mp[ 



(5) 



as well as the Klein- Gordon equations 

(j) + 3ii H — = , 

ip + SHyj H — = . 



(6) 



where H = a/a is the Hubble parameter, a is the scale 
factor and a dot denotes derivative with respect to cosmic 
time. 

In order to study the two-flelds dynamics of the hybrid 
model, without assuming slow-roll, one has to integrate 
these equations numerically from a given set of initial 
conditions (IC) for the fleld values. Throughout the pa- 
per we will deflne a successful IC as a point in fleld space 
that lead to a sufficiently long phase of inflation to solve 
the horizon and flatness problem. We will assume that 
TV = ln(a/aini) — 60 e-folds is the critical value required, 
though this value can change by a factor of two depending 
on the reheating temperature and the Hubble parameter 
at the end of inflation [H,!!!]. However, generically, once 
inflation starts it lasts for much more than 60 e-folds and 
our results are not sensitive to the peculiar value chosen. 



C. The set of successful initial field values 

As already mentioned in the introduction, the space of 
successful IC for the fleld values alone has been discussed 
in Ref. [l^ and found to be composed of a intricate en- 
semble of points organized into continuous patterns. In 
Fig. [Tl we have represented the mean number of e-folds 
generated at each sub-Planckian initial fleld values, for 



* Throughout the paper, rUpi denotes the physical Planck mass, 
and Mpi stands for the reduced Planck mass Mpi ~ 0.2mpi ~ 
2.4 X IQi** GeV. 




FIG. 1: Mean number of e-folds obtained from 512 initial 
field values in the plane (-i/Ji/Mpi, ^i/Afpi). This figure has 
been obtained by averaging the number of e-folds (truncated 
at 100) produced by 2048"^ trajectories down to 512^ pixels. 
The potential parameters have been set to M = O.OSmpi, 



■ 636mi, 



3 X 10"*m^ 



a set of fixed potential parameters and assuming vanish- 
ing initial velocities. We have computed the trajectories 
obtained from 2048^ initial fleld values and stopped the 
integration when the flelds are trapped in one of the min- 
imum of the potential, i.e. for < V/{3M^,), or when 

the accelerated expansion exceeds 10^ efolds. The re- 
sulting grid has a small intricated structure of successful 
regions spread over the whole plane which ends up being 
difficult to represent in a figure. As a result, we have 
chosen to present in Fig.[T]a downgraded 512^ pixels im- 
age in which each pixel has been given a color according 
to the average number of e-folds obtained in our original 
2048^ grid. A given pixel may therefore hide both suc- 
cessful and unsuccessful initial field values and the color 
measures their relative density. An higher resolution im- 
age would be self-similar to Fig. [1] with more thinner 
successful domains visible. 

Notice that we recover the inflationary valley as the 
white vertical narrow strip located along ijji — whereas 
the minima of the potential are along the horizontal axis 
a,t ip — ±G.15Mpi (for M = G.0377ipi as chosen in the flg- 
ure) . The black region in Fig. [T] precisely corresponds to 
the trajectories "below" the critical point 4> < (pc which 
are fast-rolling inside the minima. In analogy with the 
anamorphosis of light produced by a distorted mirror, 
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each point outside the inflationary valley is connected by 
a trajectory to a point inside the inflationary valley. The 
trajectory flrst fast-rolls towards the bottom of the po- 
tential, and after a few rebounds becomes oriented along 
the valley, climbs it and then produces inflation when 
slow-rolling back down. There is thus a one-to-one corre- 
spondence between the IC and the point in the valley for 
which the trajectory stops to climb and starts to slow- 
roU. 

It was shown in Ref. [l^l that such "anamorphosis 
points" can cover up to 15% of the total area when re- 
stricting the IC to sub-planckian values. Moreover, as 
can be checked in Fig. [1] these regions exhibits a frac- 
tal looking aspect. Before studying the influence of the 
potential parameters and initial field velocities, one may 
wonder if the area of this two-dimensional set of points is 
indeed well-defined? Equivalently, do new successful re- 
gions appear inside unsuccessful domains and conversely? 
In order to quantify how much the anamorphosis points 
are a probable way to have inflation in the whole pa- 
rameter space, we first address the question of defining 
a measure on the initial field values space. In particular, 
this requires to determine the dimension of the set 

5 = {(0i,V^i)/iV>6O}. (7) 

D. Chaotic dynamical system 

1. Phase space analysis 

As suggested by Fig. [1] at fixed potential parameter 
values, the dynamical system defined by Eqs. ^ and 
seems to exhibit a chaotic behaviour. In particular, 
the sensitivity of the trajectories to the initial field values 
comes from the presence of three attractors. Two of them 
are the global minima of the potential, A4± respectively 
a,t {(j) — Q,ip = ±M), in which all classical trajectories 
will end, whereas the less obvious is a quasi-attractor X 
defined by the inflationary valley itself (V' = 0, > 0c)- 
Indeed, whatever the initial field values, as soon as the 
system enters slow- roll one has (in Planck units) , 

where ei is the first Hubble flow function [s^. The sys- 
tem therefore spends an exponentially long amount of 
cosmic time in the valley. The sensitivity to the initial 
conditions comes from the presence of these three attrac- 
tors: either the trajectory ends rapidly into one of the 
two minima, or it lands on the valley where it freezes. 

A phase space plot is represented in Fig. [2] in which 
we have computed 25 trajectories from a grid of initial 
field values. The inflationary valley clearly appears as the 
attractor with quasi null velocity vector (ei <C 1), while 
around the two global minima, two "towers" appear due 
to the field oscillations around them. 




FIG. 2: Phase space ii^ ((/>,?/)) for 25 trajectories and vanish- 
ing initial velocities. The potential parameters are fixed to the 
values M = 0.03mpi,/i = 636mpi,i/ = 6.36 x 10"*. All tra- 
jectories end on the three attractor of the dynamical system: 
the two global minima of the potential, and the inflationary 
valley with almost vanishing slow-roll velocity. These three 
attractors induce the chaotic behavior. 



2. Basins of attraction 

From the definition of S in Eq. (O, one has 

S^F-\T), (9) 

where F((j}, ip) stands for the mapping induced by the 
differential system of equations (O and The set of 
successful initial field values S is therefore the basin of 
attraction of the attractor I [13, • Since the attractor 
X is a dense set of dimension 2 and F is continuous, one 
expects S to contain a dense set of dimension 2 [sH] . As 
can be intuitively guessed, the boundary of S can how- 
ever be of intricate structure because of the sensitivity 
to the initial conditions: two trajectories infinitely close 
initially can evolve completely differently. As we show in 
the following, S is actually a set of dimension two having 
a fractal boundary of dimension greater than one. 

Finally, by the definition of a continuous mapping, 
all parts of S, boundary included, must be connected 
together and to the inflationary valley I. The fractal 
looking aspect of Fig. [T] is only induced by the intricate 
boundary structure of S which is exploring all the ini- 
tial field values space. The fractality of the boundaries 
of the space of initial fields values was first mentioned 
in Ref. [33], but the study was restricted to a small re- 
gion of the field space and the model included dissipative 
coefficients. As an aside remark, let us notice that the 
existence of a fractal boundary may have strong impli- 
cations in the context of eternal chaotic infiation: there 
would exist inflationary solution close to any initial field 
values. 

In order to quantify the chaotic properties of the dy- 
namical system defined by the mapping V'), we turn 
to the calculation of the Lyapunov exponents. 
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FIG. 3: Highest Lyapunov exponent as a function of the ini- 
tial field values in the orginal hybrid model. The potential 
parameters are the same as in Fig. [1] The field evolution is 
therefore stable on the set 5 of successful initial field values 
(black) but exhibits chaotic behaviour elsewhere. 



3. Lyapunov exponents 



the differential system of Eqs. ([5]) and ([6]). If the largest 
exponent is positive, then the invariant set is chaotic. 

In Fig. [3l we have computed the largest Lyapunov ex- 
ponent at each point of the plane ((^j, The numerical 
method we used is based on Refs. [40, |4l[ and uses the 
public code LESNLS. Let us notice that since X is only a 
quasi- attractor, we have stopped the evolution at most 
when = F/(3Mp[), i.e. just before the fields would 
classically enter either or A^_. As can be seen, all 
points belonging to S exhibit the same and small nega- 
tive Lyapunov exponent: the invariant set S is therefore 
non-chaotic. On the other hand, all the other initial field 
values associated with the basins of attraction of A^± 
have a positive Lyapunov exponent. For those, the field 
evolution is chaotic and exhibits a sensitivity to the ini- 
tial conditions. Notice that these exponents may slightly 
vary from point to point due to our choice to stop the in- 
tegration at iJcnd instead of the classical attractors M±. 
This is particularly visible for the trajectories starting 
close to iJcnd (green shading): there is not enough evolu- 
tion to get ride of the transient evolution associated with 
the initial conditions. 



E. Fractal dimensions of S and its boundary 

1. Hausdorff and box-counting dimension 

Since we suspect a set with fractal properties, the nat- 
ural measure over S, extending the usual Lebesgue mea- 
sure, is the Hausdorff measure. The s-dimensional Haus- 
dorff measure of S is defined by [s^ 



The Lyapunov exponents at an initial point Xi — 
((/), "0, (/<Ar, '0^Ar)|i measures how fast two infinitely close 
trajectories mutually diverge or converge. They give a 
mean to characterize the stretching and contracting char- 
acteristics of sets under the mapping induced by the dif- 
ferential system. A small perturbation dx around the 
trajectory xi^) '^il^ evolve according to 



dSx 

dN 



^dF-SX: 



(10) 



where dF stands for the Jacobian of the differential sys- 
tem F. The Lyapunov exponents at the initial point 
Xi and along the direction Sxo are the numbers defined 
by M 



lim — In — — — 

N^^ N \Sxo\ 



(11) 



where Sx{N) is the solution of Eq. pop with dxiO) — 
dxo and x(0) = Xi- If the considered set is an attrac- 
tor or an invariant set of the differential system having a 
natural measure, one can show that the exponents do not 
depend on the initial point Xi- At fixed potential param- 
eters, there are four Lyapunov exponents associated with 



n%S) = lun inf <j J2 \U^\' / 5 C Q (7, ; \U,\ < S 



. i=l 



i=l 



(12) 

In this definition, the sets Ui form a 5-covering of S 
and the diameter function has been defined by \U\ = 
svip{\x — y\ / x,y € U}. As a result, 1-L'^{S) is the small- 
est sum of the sth powers of all the possible diameters 
5 of all sets covering 5, when (5^0. Having such a 
measure, the fractal dimension of S is defined to be the 
minimal value of s such that the Hausdorff measure re- 
mains null (or equivalently the maximal value of s such 
that the measure is infinite). In practice, measuring the 
Hausdorff dimension using this definition is not trivial, 
due to the necessity of exploring all (5-coverings. How- 
ever, in our case, we are interested in the fractal proper- 
ties of a basin of attraction associated with a continuous 
dynamical system and one can instead consider the so- 
called box-counting dimension [s^ . This method simply 
restricts the class of the Ui to a peculiar one, all having 
the same diameter 5. When the mapping F is self-similar, 
one can show that box-counting and Hausdorff dimen- 
sions are equal. In general, the Hausdorff dimension is 
less or equal than the box-counting one. Here, F being 
a contracting continuous flow, we expect the equality to 
also hold. 
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To define the box-counting dimension, we cover the 
set S with grids of step size S, and count the minimal 
number of boxes N{6) necessary for the covering. The 
box-counting dimension is then given by 

= Km . (13) 

<5-,0 log(l/(5) ^ ^ 

This method has the advantage to be easily implemented 
numerically and, in the following, we will apply it to cal- 
culate the dimension of S and its boundary. 

2. Fractal boundary of S 

For each randomly chosen point of the plane (i/iijT/'i), 
we compute three trajectories. The first one starts from 
the point under consideration while the two others have 
initial conditions modified by +S and —d along one direc- 
tion (for example along (f>, but the chosen direction does 
not affect the result). For each of these trajectories, we 
determine in which attractor {A4± or X) the flow ends. 
Since we are interested in the boundary of S, we calcu- 
late the proportion f{6) of points for which at least one 
trajectory ends in T, and another in A^_|_ or The 
process is iterated for increasingly smaller values of 6 and 
we evaluate how the area of the 5-grid covering of S scales 
with S. So strictly speaking, our evaluation of the box- 
counting dimension is made through the determination 
of the Minkowski dimension of the boundary of S fs^ . 
From Eq. (|13p . assuming that, at small 6, 

fiS) (X 5", (14) 

the box-counting dimension of the S boundary is then 
given by [S^l 

DB^2~a. (15) 

In Fig. 21 we have plotted f{5) as a function of S 
at fixed potential parameters. We recover the expected 
power law behaviour, the slope of which is approximately 
a ~ 0.80. As a result, the boundary of S is indeed a frac- 
tal of box-counting dimension 

Db =i 1.20. (16) 

Notice that this value depends on the chosen set of po- 
tential parameters, as one may expect since they affect 
the shape of S and the typical size of the structures. 

3. Dimension of S 

In order to determine the box-counting dimension of S 
itself one can apply a similar method than the one used 
for its boundary. Now f{S) denotes the proportion of 
points for which at least one of the three trajectories end 
in the attractor X (this condition therefore includes also 
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FIG. 4: Fraction of initial field values in a 5-sized box in- 
tercepting the set S, as a function of 5. The field has been 
restricted to sub-planckian values and the potential parame- 
ters are fixed to A = A' = 1, m = 10~®mpi and M = O.OSrripi. 
The exponent a of the power law dependency gives the box- 
counting dimension Db = 2 — a ~ 1.2 showing that 5 pos- 
sesses a fractal boundary. 



0.30 - 



0.20 - 



10-' lO-"* 10-5 10-" 0.001 0.01 

<5 

FIG. 5: Fraction of initial field values leading to inflation in a 
5-sized box as a function of S. The potential parameters are 
the same as in Fig. U) Once the box is small enough to be 
fully contained in 5, f{S) remains constant. As a result, the 
box-counting dimension of S is Db = 2 and the interior of S 
is not fractal. 



the points belonging to the boundaries). The resulting 
power-law is represented in Fig. [5] 

For small enough values of S, the (S-sized boxes are 
small enough to be fully contained in S and the function 
f{5) appears to be constant in that case. As a result, 
the box-counting dimension of 5 is 2. We therefore con- 
clude that, like for the well-known Mandelbrot set [43 |. 
the boundary of S is fractal but the set of successful infla- 
tionary points is not and has the dimension of a surface. 
Consequently, although the boundary of S has an infinite 
length (Db = 1-2), it has a vanishing area: the Hausdorff 
dimension of S (boundary included) is therefore also 2. 
As a result, the two-dimensional Hausdorff measure on 
S reduces to the usual two-dimensional Lebesgue mea- 
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sure and this will be our choice for defining a probability 
measure in the rest of the paper. 

As previously emphasized, the potential parameters 
and initial field velocities have been fixed in this section 
and the set S is actually the two-dimensional section of 
an higher dimensional set, whose boundary is also cer- 
tainly fractal (and therefore of null measure). Since one 
can no longer use griding method to explore such a high 
dimensional space, we move on in the next section to a 
MCMC exploration of the full parameter space to assess 
the overall probability of getting inflation in the hybrid 
model. 



MCMC elements directly sample the posterior probabil- 
ity distribution Tr{6\d) of the model given the data. 

In our case, we can similarly define a likelihood £ as a 
binary function of the potential parameters, initial field 
values and velocities. Either the trajectory ends up on X 
and produces more than 60 e-folds of inflation, or it does 
not. In the former case we set £ = 1 whereas £ = for 
no inflation. The function tt we sample is then defined 
by TT = CP{9) where 9 stands for field values, velocities 
and potential parameters and P is our prior probability 
distribution that we will discuss in the next section. 



III. PROBABILITY DISTRIBUTIONS IN 
HYBRID INFLATION 

The aim of this section is to use Monte-Carlo-Markov- 
Chains (MCMC) techniques in order to explore the whole 
parameter space, including the initial field velocities and 
all the potential parameters. With unlimited computing 
resources, we could have used a griding method to lo- 
calise the hypervolumes in which inflation occurs, as we 
have done for the two-dimensional plane in the 

previous section. For the original hybrid model, we have 
in total seven parameters that determine a unique tra- 
jectory: two initial field values, initial field velocities and 
the three potential parameters M, /i and v. To probe 
this seven-dimensional space, more than just measuring 
the hypervolume of the successful inflationary regions, 
we define a probability measure over the full parameter 
space. Using Bayesian inference, one can assess the pos- 
terior probability distribution of all the parameters to 
get enough e-folds of inflation. Monte-Carlo-Markov- 
Chains (MCMC) method is a widespread technique to 
estimate these probabilities, its main power being that 
it numerically scales linearly with the number of dimen- 
sions, instead of exponentially. 

Several algorithms exist in order to construct the 
points of a Markov chain, the Metropolis-Hastings al- 
gorithm being probably the simplest 0, [3] . Each point 
Xi+i, obtained from a Gaussian random distribution (the 
so-called proposal density) around the previous point Xi , 
is accepted to be the next element of the chain with the 
probability 



Pix^+l) 



1, 



7r(a;i+i) 



(17) 



where Tr(x) is the function that has to be sampled via the 
Markov chain. MCMC methods have been intensively 
used in the context of CMB data analysis [H, [H, \^ HI, 
Eil where the function 7r(e'|d) cx C{d\e)P{e) is the pos- 
terior probability distribution of the model parameters 
given the data. In the context of Bayesian inference, this 
one is evaluated from the prior distributions P{9) and the 
likelihood of the experiment C{d\9). After a relaxation 
period, one can show that Eq. (fT7|) ensures that tt is the 
asymptotic stationary distribution of the chain [soj . The 



A. Prior choices 

MCMC methods require a prior assumption on the 
probability distributions of the fields, velocities and po- 
tential parameters. As we only consider in this work 
the initial conditions and parameters space leading to at 
least 60 e-folds of inflation, the prior choices are only 
based on theoretical arguments. These arguments can 
be linked to the framework from which the potential is 
deducted. If one considers the hybrid model to be em- 
bedded in supergravity, the fields have to be restricted to 
values less than the reduced Planck mass. We adopt here 
this restriction for initial field values, not only because 
of this argument, but also because it has been shown in 
Ref. j22,] that if super-plankian fields are allowed, trajec- 
tories become generically successful. On the other hand, 
the model was considered to suffer some fine-tuning when 
one of the fields has to be order of magnitudes smaller 
than the other. As inflation is not possible for very small 
initial values of both fields (because of the Higgs instabil- 
ity) , we have considered a flat prior for initial field values 
in [— Afpi, Mpi] as opposed to a flat prior for the logarithm 
of the flelds. Note that one has to include negative values 
of the fields in order to take into account the orientation 
of the initial velocity vector. 

Concerning the initial field velocities, from the equa- 
tions of motion, one can easily show that there exist a 
natural limit^ 



dN 



diA 
dN 



< 6. 



(18) 



Similarly, our prior choices are flat distributions inside 
such a circle in the plane {(j),N, '4>,n), where ", TV" denotes 
a partial derivative with respect to the number of e-folds. 

In the absence of a precise theoretical setup determin- 
ing the potential there are no a priori theoretical con- 
straint on its parameters M , /i and v. Let us just mention 
that for /i < 0.3, the dynamics of inflation in the valley 
is possibly strongly affected by slow-roll violations '2^1. 
As a result, with the concern to not support a particular 



This is just the limit ei < 3 in Planck units [35 
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mass scale, we have chosen the fohowing flat priors on 
the logarithm of the parameters: 

-1 < log-^ < 4, 
mpi 

-3<logl^-<-0.7, (19) 

-6 < log < 2, 

in which the upper and lower limits have been set for 
numerical convenience, and M < Mpj. 

Notice that the A dependencies are not important here 
because this parameter only rescales the potential and 
thus does not change the dynamics. 

In the next sections, we perform the MCMC ex- 
ploration of the parameter space from these priors. 
Firstly by reproducing the results of Sec. in the two- 
dimensional section then by including the initial 
field velocities and finally by considering all the model 
parameters. Unless otherwise mentioned, the chains con- 
tain 10® points, which corresponds to one percent error 
on the marginalised probability distributions. In the fig- 
ures, the overall values of the posterior probability den- 
sity distributions have not been represented since they 
are determined by the imposing the integral over the pa- 
rameters to be equal to one. 




-0.2 -0.1 0.1 0.2 
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FIG. 6: Two-dimensional posterior probability distribution in 
the plane {<j>i,'4>i) leading to more than 60 e-fold of inflation 
in the hybrid model. Notice that its integral over the plane 
is normalised to unity. The dark blue regions corresponds to 
a maximal probability density whereas it vanishes elsewhere. 
The potential parameters are set to M = 0.03 rripi, = 6.36x 
10~^ rripi, fi = 636 nipi. As expected, the MCMC exploration 
matches with the griding methods (see Fig. [TJ . 



B. MCMC on initial field values 

In order to test our MCMC, we have first explored the 
space of initial field values leading or not to more than 60 
e-folds of inflation. The potential parameters have been 
fixed to various values already explored by griding meth- 
ods in Sec. HIl and Ref. [l^l, while the initial velocities 
are still assumed to vanish. The MCMC chain samples 
have been plotted in Fig. [6l Notice that to recover the 
fractal structure of the boundary of S, one has to adjust 
the choice of the Gaussian widths of the proposal density 
distribution. If those are too large, the acceptance rate 
will be small because the algorithm tends to test points 
far away from the last successful point, and if they are 
too small the chains remain stuck in the fractal struc- 
tures without exploring the entire space. Nevertheless, 
with an intermediate choice. Fig. [S] shows that the intri- 
cate structure of the boundary of S can be probed with 
the MCMC. More than being just an efficient exploration 
method compared to griding, the MCMC also provides 
the marginalised probability distributions of and ^pi 
such that one gets inflation. They have been plotted 
in Fig. [in] (top two plots), the normalisation being such 
that their integral is unity. As one can guess from Fig. [51 
with vanishing initial velocities and a fixed set of poten- 
tial parameters, infiation starting in the valley is not the 
preferred case since the area under the distribution of 
-01 outside of the valley is larger than inside. Moreover, 
these distributions take non- vanishing values everywhere 



and there is therefore no fine-tuning problem. Of course, 
one still have to consider the other parameters and this 
is the topic of the next sections. 



C. MCMC on initial field values and velocities 

The initial values of the field velocity are inside a 
disk of radius -\/6 in the plane {4>^n,4',n) (in Planck 
units). The marginalised two-dimensional posteriors for 
the initial field values is plotted in Fig. [7| whereas the 
marginalised posterior for each field are represented in 
Fig. [TUj (middle line). Even if non- vanishing velocities 
are considered, the successful infiationary patterns re- 
main. Notice that they appear to be blurred simply be- 
cause of the weighting induced by marginalising the full 
probability distribution over the initial velocities. 

In Fig. [51 we have also represented the marginalised 
posterior probability distribution for the modulus and 
direction of the initial velocity vector. Their flatness im- 
plies that there are no preferred values. This is an im- 
portant result because one could think that large initial 
velocities could provide a way to kick trajectories in or 
out of the successful regions. This actually never happens 
because of the Hubble damping term in the Friedmann 
equations, allowing only a generation of a small number 
of e-folds before the trajectory falls in one of the three 
attractors. 
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FIG. 7: Two-dimensional marginalised posterior probability 
distribution for the initial fields values. The marginalisation 
is over the initial field velocities whereas the potential pa- 
rameters are still fixed. The shading is proportional to the 
probability density value while the two-dimensional integral 
over the plane is equal to one. Although the inflationary valley 
has the highest probability density, its area remains restricted 
such that the most probable initial field values to get inflation 
are still out of the valley (see Fig. [TOj . 



D. MCMC on initial field values, velocities and 
potential parameters 

The most interesting part of the exploration by MCMC 
technique concerns the study of the full parameter space. 
The only restriction being associated to the necessity of 
M < Mpi as discussed in Sec. IIII Al The chains contain 
200000 elements and the estimated error on the posteriors 
is about a few percents. 

We have plotted in Fig. [9] the marginalised two- 
dimensional posterior for the initial field values. In com- 
parison with Fig. [5] and [71 the most probable initial field 
values are now widespread all over the accessible values; 
the intricate patterns that were associated with the suc- 
cessful field values (at fixed potential parameters) are 
now diluted over the full parameter space. The resulting 
one-dimensional probability distributions for each field 
are plotted in Fig. [10] (bottom panels) . One can observe 
that the ?/; distribution is nearly flat outside the valley but 
remains peaked around a extremely small region around 
■0 = 0. Integrating over the field values, initial conditions 
outside the valley are still the preferred case. 

Concerning the probability distributions of the modu- 
lus V and the angular direction 9 of the initial velocity 
vector, results integrated over the whole parameter space 
do not present qualitative differences compared to the 
posteriors with fixed potential parameters, as one may 
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FIG. 8: Marginalised posterior probability distributions for 
the modulus (top) and angle (bottom) of initial field velocity. 
The thin superimposed blue (lighter) curves are obtained at 
fixed potential parameters, while the thick red are after a full 
marginalisation over all the model parameters. As expected 
from Hubble damping, all values are equiprobable since the 
field do not keep memory of the initial velocity. 

expect since the Hubble damping prevents the initial ve- 
locities to influence the dynamics (see Fig. [5]) . 

The marginalised probability distributions for the po- 
tential parameters are represented in Fig. [11] These pos- 
teriors seem to indicate that the three parameters are 
bounded but one should pay attention to the influence 
of our prior choices over the posterior jsij . In fact, the 
posteriors for M and v are found to depend on our prior 
choice: changing the upper or lower limit on the i/-prior 
(or M-prior) affects the values at which the M- and v- 
posteriors fall off. Such a situation is typical of the exis- 
tence of correlations between these two parameters. We 
have therefore computed the two-dimensional posterior 
distribution in the plane {v^ , M) and found out that this 
probability distribution clearly exhibits a correlation be- 
tween these two parameters: the lower bound on M de- 
pends on the minimal allowed value of v in the prior. 
Such a correlation comes from the fact that, to realize 
enough inflation for a given IC, the critical/instability 
point 0c should not be larger than the anamorphosis im- 
age of the IC (the point in the valley where slow-roll 
starts). Restricting initial fields to sub-Planckian values 
leads to an uper bound on the largest image, and thus an 
upper bound on the instability point. From Eq. ((H]), the 
relevant quantity that is constrained is the combination 
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occur. As discussed in the introduction, the original hy- 
brid model under scrutiny is however a toy model known 
to be disfavored by the current CMB data. In this re- 
spect, one may wonder whether our results are peculiar 
to this model or can be generalised to other more realis- 
tic two field inflationary models. This point is addressed 
in the next section in which we have performed a com- 
plete study of the SUGRA F-term hybrid inflation. In 
that model, the dynamics depends on only one potential 
parameter; also constrained by cosmic strings formation. 
The challenge will thus be to confront this constraint to 
the natural bounds that can be deducted from MCMC 
methods by requiring enough e-folds of inflation. 



IV. PROBABILITY DISTRIBUTIONS IN 
F-SUGRA INFLATION 



»"pi 



FIG. 9: Two-dimensional marginalised posterior probability 
distribution for the initial fields values. The marginalisation 
is over the initial field velocities and all the potential param- 
eters. The shading is proportional to the probability density 
value. The infiationary valley is still visible around ^pi — 
and the posterior takes non-vanishing values everywhere in 
the {4)i,^i) plane. 



The minimal supersymmetric versions of hybrid infla- 
tion are known as the F-term and D-term inflationary 
models [H, [13, [III , where the slope of the valley is gener- 
ated by radiative corrections. The F-term model is com- 
patible with the current CMB data since a red spectrum 
of the cosmological perturbations is generic 
In addition, this model is more predictive and testable 
than its non-SUSY version since it contains only one cou- 
pling constant and one mass scale. 



We have plotted in Fig. [12] the marginal posterior dis- 
tribution associated with the parameter \og{'\/2v'^ / M) , 
and at 95% of confidence level, we find 

< 4 X 10-3. 20 

M ^ ' 

The parameter is the other constraint that the 
MCMC exibits. It is explained by the possible appari- 
tion of slow-roll violations in the valley, when ji becomes 
too small. These slow-roll violations prevent the genera- 
tion of an inflationary phase if the trajectory climbs too 
high in the valley. At a two-sigma level, one has 

-^>1.7. (21) 

This lower limit is equivalent to the upper limit on m 
observed in : a large inflaton mass induces a fast roll 
evolution and requires super-planckian initial conditions 
to realize inflation in a chaotic way. Let us stress that 
these constraints come only from requiring enough infla- 
tion in the hybrid model whatever the initial field values, 
velocities, and other potential parameters. In this re- 
spect, the limits of Eqs. ([^0]) and (PT|) can be considered 
as "natural" . 

To conclude this section, we have shown that inflation 
is generic in the context of the hybrid model and we have 
derived the marginalised posterior probability distribu- 
tions of all the parameters such that 60 e-folds of inflation 



A. The model 

In the following, we are analysing the space of initial 
conditions and model parameters leading to enough in- 
flation for the so-called F-term model based on the su- 
perpotential [l3l 

<fl = «:5(<i>+$_ - Af2) . (22) 

The inflaton is contained in the superfield S. The Higgs 
pair <I>4-,<I'_ is charged under a gauge group G, that is 
broken at the end of inflation when the Higgs pair de- 
velop a non-vanishing expectation value (vev) M . The 
superpotential leads in global SUSY to a tree level po- 
tential 

where the effective inflaton s and Higgs fleld can be 
made real and canonically normalized [s = \/2 3?('5'), 4> — 
25R((/)+) = 23fi(0„)]. The local minima of the potential at 
large S provide a flat direction for the inflaton s: Vq — 

This tree level flat direction is lifted by two effects. 
Firstly, radiative corrections are induced by the SUSY 
breaking that supports inflation. In addition, if the 
field values are close to the reduced Planck mass Mpi, 
one should expect supergravity corrections S/Mpi to the 
tree level potential. The radiative corrections along the 
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FIG. 10: Marginalised posterior probability distributions for the initial field values <j)i and il'i- The top panels correspond to 
vanishing initial velocities and fixed potential parameters, the middle ones are marginalised over velocities at fixed potential 
parameters, while the lower panels are marginalised over velocities and all the potential parameters. 



inflationary valley can be derived using the Coleman- 
Weinberg formula 54]. They reduce to [I4I 



v-rioop(s) 



327r2 



21n- 



A2 



+ (z + l)2ln(l + z-i) 



+ {z-lf\n{l- z-^) 



(24) 



that effectively^ A/" = 2, 3. For the sake of generality, we 
will assume that M can take values in the range [2, 126]. 
This model is also known to generically produce cosmic 
strings at the end of inflation [2^ and this irnposes an 
upper limit on the inflationary mass scale [ssl [Sa . [s^ 

126 

Af < 2 X 15 GeV, k < 7 x lO"'^ -jj- . (25) 



where z = /M"^^ M stands for the dimensionality of 
the representations to which <f>± belong and A represents 
a non-physical energy scale of renormalization. Realistic 
values of M can be derived from the embedding of the 
model in realistic SUSY Grand Unified Theories (GUT) 
as shown in Ref. [2^. For example, in the case of an 
embedding of the model in SUSY SO(IO), <i>± belong 
to the representation 16, 16 or 126, 126. However, as 
pointed out in Ref. [SSI ]. it is possible that only some 
components of take a mass correction of order M so 



Secondly, SUGRA corrections also contribute to lifting 
the tree-level flat direction and will be taken into account 
since the field values we are probing are not always negli- 
gible compared to the Planck mass. It has been noticed in 
Ref. [ll] that the F-term hybrid inflation model doesn't 
suffer from the ry-problem only when the Kahler potential 



" This depends on the mass spectrum of the assumed GUT model. 
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FIG. 11: Marginalised probability distribution for the poten- 
tial parameters of the hybrid model. Notice that some of the 
bound are set by the prior choices. 



is (close to) minimal 



if ~ |5|2 + |$+|2-f|$_|2 , 



(26) 



which is what we assume in the following. In terms of the 
canonically normalized effective inflaton s and waterfall 



We will restrict ourselves to minimal SUGRA corrections, ne- 
glecting SUSY breaking soft terms and the non-renormalizable 
corrections to the superpotential (see [s^ . Issl . l55j for an analysis 
of their effects). 




FIG. 12: Prior independent marginalised posterior probability 
distribution for the parameter u^/{MMp-i). This parameter 
fixes the position of the instability point and a too large value 
may prevent infiation from occurring in the sub-planckian 
field regime [see Eq. [20)) ] . 



fields ^p, the SUGRA corrected potential now reads 




1 - 



2M2, 



(27) 



The dynamics along the inflationary valley is driven by 
the radiative corrections and by the SUGRA corrections. 
The radiative corrections play a major role in the last e- 
folds of inflation (thereby generating the observable spec- 
tral index), whereas most of the dynamics takes actually 
place for fleld values dominated by the SUGRA correc- 
tions. We have calculated the amplitudes for both correc- 
tions and found that only at the end of the inflationary 
potential (for s e [M, 8M] if TV = 3 and s G [M, 3.5A/] if 
TV = 126), the radiative corrections may dominate over 
the SUGRA corrections. In the present work, the regions 
of the parameter space leading to inflation do not depend 
on the very last part of the field evolution: as soon as 
60 e-folds are obtained, the initial conditions are consid- 
ered successful and this generically occurs in the valley 
at larger fleld values. Outside the inflationary valley, we 
therefore expect the tree level dynamics to dominate over 
the radiative corrections, especially for small coupling k. 
There also, in addition to the tree level, at large fields, 
SUGRA corrections are expected to be important. 

Resulting of these considerations, we have neglected 
radiative corrections and used for our study below the 
potential of Eq. (P7)) . 
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Values of M 


Area of 5 (%) 


M = 10~^ mpi 
M = 10"^ mpi 
M = 10"^mpi 
M = IQ-^mpi 


(exact) 
12.9 ± 0.1 
12.0 ± 0.3 
10.3 ±0.5 



TABLE I: Percentage of successful initial field values, at van- 
ishing initial velocities, for various values of the potential pa- 
rameter M. The error bars come from the finite numerical 
precision, which decreases with M. 
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FIG. 13: Mean number of e-folds (truncated at 100) obtained 
from 512^ initial field values (-^i/Afpi, si/Mpi) for the SUGRA 
F-term model. The initial field velocities are assumed to van- 
ish and the potential parameter is fixed at M = 10~^mpi. 
As for the original hybrid model, the mean is computed from 
2048^ trajectories (see Fig. [T]). The set of initial field values 
producing enough inflation is again of dimension two with a 
fractal boundary. 



log (5 



FIG. 14: Fraction of initial field values in a 5-sized box 
intercepting the set 5 as a function of 5 for the SUGRA 
F-term model. The potential parameter has been fixed to 
M — 10~^mpi. The box-counting dimension of boundary of 
S is given by the power law behaviour for small 5 and found 
to be Db — 1.5. 



B. Fractal initial field values 

The analysis of the SUGRA F-term model of inflation 
has been conducted along the lines described in Sec. HIl 
and Sec. mil We have flrst verified that, at fixed potential 
parameter M and vanishing initial velocities, the set of 
initial field values S defined by Eq. ([7]) is two-dimensional 
with a fractal boundary. In Fig. 1131 we have represented 
the set S of successful initial field values for the mass 
scale M = 10~^mpi. Notice that the coupling constant 
K, being an overall factor, it doesn't impact the dynamics 
of the fields. Our study is therefore valid for any value 
of K, and of the dimensionality of the Higgs field A/", since 
the relationship M{k) depends only on M . 

As for the original hybrid model, the highest Lyapunov 
exponent for the successful initial field values is negative 
and the set S is non-chaotic. Outside of 5, trajectories 
have positive Lyapunov exponents and exhibit chaos. 

For vanishing initial velocities, we have reported in Ta- 
bleUthe area occupied by the set S in the plane (si,V'i) 
for various section along the potential parameter M . Like 
for the original hybrid model, we recover a significant pro- 
portion of successful initial field values outside the valley. 
This result holds even for M ^ 1 though at small M, the 



potential becomes very flat and the number of oscilla- 
tions of the system before being trapped in the inflation- 
ary valley can exceed 10^. Simulations become therefore 
more time-consuming and error-bars in Tab. U increase. 
Reducing M also reduces the typical size of structures in 
the plane (si,?/'i), which evolves from Fig. [T3]to a more 
intricated space of thinner successfuU IC. As suggested 
by the Tab. U we will see below that this doesn't affect 
the probability of getting inflation by starting the field 
evolution outside the valley. 

Concerning the fractal properties of 5, we have ap- 
plied the same method as in Sec. HIE II to compute the 
box-counting dimensions of S and its boundary. As ex- 
pected, we recover that S is of box-counting dimension 
two whereas the function f{5) for its boundary is repre- 
sented in Fig. [Ml We obtain that, as in the non-SUSY 
case, the boundaries are fractal with dimension 

Db ~ 1.5 . (28) 

These results allow us to use the usual Lebesgue mea- 
sure to define the probability distribution over the whole 
parameter space. 
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FIG. 15: Marginalised posterior probability distributions for 
the initial field values (upper panels) and the initial velocities, 
modulus V and angle 9. The F-SUGRA inflationary valley 
has a slightly higher probability density around = but is 
extremely localised: as a result, inflation is more probable by 
starting out of the valley. 



C. MCMC on the initial field values, velocities and 
the potential parameter 

As already mentioned, there is only one potential pa- 
rameter M in F-term SUGRA model that may influence 
the two field dynamics. The goal of this section is to 
evaluate the probability distributions of the initial field 
values, velocities and of M such that inflation lasts more 
than 60 e- folds. As for the original hybrid model, we have 
performed an MCMC analysis on the five-dimensional 
parameter space defined by Sj, w, and M where 



ds 
dN 



— V cos 6*, 



dN 



V smf 



(29) 



We have chosen the same sub-planckian priors for the 
initial field values and initial velocities than in Sec. IIIII 
Since the order of magnitude of M is not known, we have 
chosen a flat prior on 



M 

2<log— <0. 



(30) 



The lower limit on M is motivated by computational 
rather than physical considerations. The resulting 
marginalised posterior probability distributions for each 
of the parameters are represented in Fig. 1151 The chains 
contain 400000 samples producing an estimated error on 
the posteriors around a few percent (from the variance 
of the mean values between different chains) . 

The posteriors for the field velocities are flat show- 
ing that all values are equiprobable to produce inflation. 
The initial field values are also flat, up to a sharp peak 
of higher probability density around ■0 = correspond- 
ing to the inflationary valley. As for the hybrid model 
of Secini after integration of these curves over the field 



FIG. 16: Marginalised posterior probability distribution of 
the mass scale M of F-SUGRA inflation. 



values, inflation is clearly more probable by starting out 
of the valley. Finally, only the posterior probability dis- 
tribution of logM is strongly suppressed at large values. 
We find, at 95% of confidence level (see Fig. [TO)). 



log(Af) < -1.33. 



(31) 



As for the original hybrid model, this limit comes from 
the condition of existence of a sub-planckian inflationary 
valley which is related to the position of the instability 
point. Indeed, from Eq. (|27|) . one finds 



dV- 



±- 



M 



"M, 



pi 



\ 



1- Wl-4 



(32) 

where we have kept only the sub-planckian solutions. 
This expression shows that there is an infiationary valley 
at xj; — Q only for M/Mpi < 1/^2, and for field values 
such that s > Sc- As a result of the two-field dynam- 
ics, we find that a valley supporting at least 60 e-folds 
of inflation require the more stringent bound of Eq. (I5T|) . 
Let us flnally notice that the most probable values we 
obtain on M to get inflation in Eq. (pij) are compati- 
ble with the existing upper bound coming from cosmic 
strings constraint: M < 10~^r7ipi (see Ref. (ssl. [56j). 



V. CONCLUSION 

In this paper, by numerically solving the two-field dy- 
namics of the original hybrid model and its SUGRA F- 
term version, we have shown that 60 e-folds of infiation 
is a generic outcome. Contrary to what is usually as- 
sumed, one does not need to fine-tune the initial field 
values around '0 = to get inflation. In fact, the infla- 
tionary valley, indeed of small extension in field space, is 
one of the three dynamical attractors of the differential 
system given by the Einstein and Klein-Gordon equa- 
tions in a FLRW universe (the others being the minima 
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of the potential). As a result, any trajectory will end 
in one of these three attractors and the set S of suc- 
cessful initial conditions therefore belongs to the basin 
of attraction of the inflationary valley. We have shown 
that such a set is connected and of dimension two while 
exhibiting a fractal boundary of dimension greater than 
one. Moreover, it occupies a significant fraction of the 
sub-planckian field regime. In order to quantify what are 
the natural field and parameter values to get inflation for 
both of these models, we have introduced a probability 
measure and performed a MCMC exploration of the full 
parameter space. It appears that the inflationary out- 
come is independent of the initial field velocities, is more 
probable when starting out of the inflationary valley, and 
favours some "natural" ranges for the potential param- 
eter values that cover many order of magnitudes. The 
only constraints being that the inflationary valley should 
at least exist. 

Let us notice that the posterior probability distribu- 
tions we have derived are not sensitive on the fractal 
property of the boundary of S. This is expected since, 
even fractal, the boundary remains of null measure com- 
pared to S. However, its existence may have implica- 
tions in the context of chaotic eternal inflation [53, . 
Indeed, the boundary itself leads to inflation and spawn 



the whole field space such that its mere existence im- 
plies that inflationary bubbles starting from almost all 
sub-planckian field values should be produced. Here, we 
have been focused to the classical evolution only and our 
prior probability distributions have been motivated by 
theoretical prejudice (flat sub-planckian prior). In the 
context of chaotic eternal inflation, our results are how- 
ever still applicable provided one uses the adequate prior 
probabilities which are the outcome of the super-Hubble 
chaotic structure of the universe [5^ . Provided the eter- 
nal scenario does not correlate with the classical dynam- 
ics, one should simply factorise the new priors with the 
posteriors presented here to obtain the relevant posterior 
probability distributions in this context. 
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